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We explore a classical instability of spacetimes of dimension D > 4. Firstly, 
we consider static solutions: generalised black holes and brane world metrics. 
Oh, The dangerous mode is a tensor mode on an Einstein base manifold of dimension 

' D — 2. A criterion for instability is found for the generalised Schwarzschild, AdS- 



Schwarzschild and topological black hole spacetimes in terms of the Lichnerowicz 
spectrum on the base manifold. Secondly, we consider perturbations in time- 



X 
U ■ 

dependent solutions: Generalised dS and AdS. Thirdly we show that, subject 



to the usual limitations of a linear analysis, any Ricci flat spacetime may be 
stabilised by embedding into a higher dimensional spacetime with cosmological 
constant. We apply our results to pure AdS black strings. Finally, we study the 
stability of higher dimensional "bubbles of nothing" . 
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1 Introduction 



Over the last few years solutions of the Einstein equations in higher dimensions have 
come to play an important role as background metrics in various physical applications. 
These range from theories of TeV gravity, where higher dimensional black holes are 
predicted to be produced in the next generation of colliders [jl], , to the gravity-gauge 
theory correspondence ||, ||] . 

Clearly the stability of such spacetimes is an important issue. One feature of 
higher dimensional spacetimes is that they often satisfy boundary conditions which 
differ from those encountered in four spacetime dimensions. This is because the 
two-sphere and two dimensional hyperbolic space are, up to discrete quotients, the 
unique Einstein manifolds in two dimensions with positive and negative curvature 
respectively. In higher dimensions there are more possibilities & [|]. These include 
metrics such as the Bohm metrics [0] that exist on manifolds that are topologically 
S d . 

In particular, we will consider the case in which the higher dimensional spacetime 
includes a d-dimensional Einstein manifold, {B,g}, which we call the base manifold, 
in a common way. In these cases we shall show that part and sometimes all of the 
stability problem may be reduced to the solution of an ordinary differential equation 
of Schrodinger form. The modes we concentrate on are transverse tracefree tensor 
harmonics on the base manifold, {B,g}. The differential equation determining sta- 
bility of the spacetime then depends on the spectrum of the Lichnerowicz operator 
on transverse traceless symmetric tensor fields of the manifold B. These modes do 
not exist in the stability analysis of, for example, the Schwarzschild black hole in four 
dimensions |, g, [ji], 0] because there are no suitable tensor harmonics on S 2 [p^ ]. 
Thus, the instabilities we discuss are inherently higher dimensional. 

Typically, the metric g on B will be such that 

R a p = e(d-l)g a p, (1) 

with e = ±1 or e = 0. This is the normalisation of S d , for example. Tildes denote 
tensors on B. 

Some examples of such spacetimes are now in order. 
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(a) Static solutions 

The spacetime is D = d + 2 dimensional and of the form 



ds 2 = -f{r)dt 2 + ^- )+ r 2 ds 2 (2) 



with 

d-1 



f(r) = e-(^j ~ -cr 2 , (3) 

and ds\ is the metric on B. The cosmological constant in d + 2 dimensions is 

R ah = c{d + l)^. (4) 

Consider first vanishing cosmological constant, c = and with e = 1 in (ffl). 
When B = S d , these are the Schwarzschld-Tangherlini black holes [|13| which are 
spatially Asymptotically Euclidean (AE). If B ^ S d one obtains generalised higher 
dimensional black holes P] which are spatially Asymptotically Conical (AC). These 
are of course not possible in four dimensions because S 2 is the only positive curvature 
Einstein manifold in two dimensions. If B = S /V where V C SO(d + 1) is discrete, 
then the spatial metric will be Asymptotically Locally Euclidean (ALE). The (d + 1)- 
dimensional Riemanian manifold with metric 

dp 2 + p 2 ds 2 d , (5) 

is called the cone C(B) with base B. It is Ricci flat precisely if g is Einstein with the 
Einstein constant of (jl])- If B ^ S d there will typically be a singularity at the vertex 
of the cone, but in our case this will often be hidden inside an event horizon. 

In forming time dependent solutions below we will replace p 2 by sin 2 p or sinh 2 p 
in (||) to obtain (d + 1) dimensional Einstein metrics of positive or negative scalar 
curvature respectively. This will have two or one singular vertices respectively. In 
the latter case the Riemannian metric will be Asymptotically Hyperbolic (AH). An- 
other way to get an Asymptotically Hyperbolic Einstein metric with negative scalar 
curvature is to take e = and replace p 2 by e 2p . 

If c is positive and e = 1, one gets a generalized Schwarzschild-Tangherlini-de 
Sitter spacetime. The static region between a cosmological event horizon and a black 
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hole event horizon is non-singular. If c is negative and e = 1, one has generalised 
Schwarzschild-Tangherlini-Anti-de Sitter, without a cosmological horizon. 

Another interesting possibility with no cosmological horizon is to take c negative 



and e = — 1 [14, 15, 16]. Now, if a = 0, the resulting metric will be singularity free. If 
B is a hyperbolic manifold B = H d /F, with T C SO(d, 1) a suitable discrete group, 
then we have an identification of Anti-de Sitter spacetime sometimes thought of as 
a topological black hole. However if B is not a hyperbolic manifold then one gets a 
singularity free topological black hole which is not locally isometric to Anti-de Sitter 
spacetime. 

Also of this form are solutions with negative cosmological constant of a sort which 



arise in brane world scenarios [17, 18|. The metric is most familiar in the form 



ds 2 = [dz 2 - dt 2 + dsfj , (6) 
with e = 0. If B is flat we obtain d + 2 dimensional Anti-de Sitter spacetime. 



(b) Time-dependent solutions 

By suitably re-interpreting our formulae we can also discuss the stability of some 
time-dependent solutions. For example a generalized D = d+1 dimensional de Sitter 
spacetime is given by 

ds 2 = -dt 2 + C ° S ^ 2 U ds 2 d . (7) 

with e = 1. This is singularity free. Changing cosh Lt to sin Lt in (0) and letting 
e = — 1 will give a generalized Anti-de Sitter spacetime which will have Big Bang and 
Big Crunch singularities at t = and t = ir/L respectively unless B is the hyperbolic 
metric on H d . 



(c) Ricci flat Lorentzian base and double analytic continuation 

A simple generalization of this time dependent situation arises if we take B to be 
a d dimensional Lorentzian Ricci flat manifold whose stability properties are known. 

For example we could consider the d + 1 dimensional Einstein manifold with 
negative scalar curvature whose metric is 

ds 2 = ^{dz 2 +dsfj, (8) 
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with e = 0, because the base is Ricci flat. If B is the flat Minkowski metric then we 
have the metric of d + 1-dimensional Anti-de Sitter spacetime. If B is a black hole 
metric, then we have black strings in Anti-de Sitter spacetime. 

Another situation in which a Lorentzian base arises is in double analytic contin- 
uation of black hole metrics. The resulting solutions describe expanding "bubbles 
of nothing" . Double analytic continuation of generalised Schwarzschild gives the 
metric 

d-l 



ds 2 



1-1- 



dr 2 

d ^ 2 + 717^1+ r2dS l (9) 



where ift is periodic and ds d is a Lorentzian metric obtained via analytic continuation 
of a Euclidean Einstein metric with e = 1. If the Euclidean base is S , then the 
corresponding Lorentzian base is just de Sitter space, dSj. 

In §2 we relate the Lichnerowicz operator on certain modes in the spacetime with 
the Lichnerowicz operator on the base manifold. This will give us equations for 
the perturbative modes. In §3 we study the stability of generalised static metrics 
by setting up a Sturm-Liouville problem. In §4 we look at perturbations in time 
dependent scenarios. In §5 we recall the Lichnerowicz spectra on some manifolds 
that give explicit examples for the results of section 3. Finally, §6 contains the 
conclusions. 

2 Lichnerowicz operator on a class of spacetimes 

2.1 A Lichnerowicz mode 

Consider a D dimensional spacetime with metric 

ds 2 D = ~f(r)dt 2 + g(r)dr 2 + r 2 ds 2 d , (10) 

where ds d is a Riemannian metric on a d = D — 2 dimensional manifold B. The 
spacetime is taken to be Einstein. 

The Lichnerowicz operator acting on a symmetric second rank tensor h is 

(A L h) ab = 2R c abd h d c + R ca h c b + R cb h c a - V c V c h ab . (11) 
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For transverse tracefree perturbations this gives the first order change in the Ricci 
tensor under a small perturbation to the metric 

9ab -> 9ab + h ab , such that h a a = V a h ab = 

Rab^ Rab + \{^ L h)ab- (12) 

The Lichnerowicz operator is compatible with the transverse, tracefree condition pC|]. 

We wish to study the stability of metrics of the form (pi]) under certain metric 
perturbations. It will be useful to have an expression for the Lichnerowicz operator 
on the spacetime in terms of the Lichnerowicz operator on the base manifold B. We 
shall impose the conditions 

ho a = h la = 0, (13) 

where 0, 1 are the t, r coordinates. These conditions are of course not a gauge choice 
and mean that we are restricting the modes we are looking at. More precisely, we are 
restricting attention to tensor modes on the base manifold and we are not considering 



scalar and vector modes. However, following [21] we argue in appendix A that, at 
least for the manifolds in which the base is compact and Riemannian with e = 1, 
the stability of the spacetime under scalar and vector perturbations is insensitive to 
the base manifold. Therefore, for these modes one may consider the base to be the 
sphere, S d . But this leaves us with just the standard Schwarzschild-Tangherlini(-AdS) 
spacetimes. These standard higher dimensional black holes are expected, although 
to our knowledge not proven, to be stable against vector and scalar perturbations. 
Therefore, we expect that it is only the tensor modes which probe the base manifold 
sufficiently to produce instabilities. Nonetheless, it would be nice to see this from an 
explicit perturbation analysis. The conditions (13) and the form of the metric ( |i~0| ) 
imply that the transeverse tracefree property of h ab (|i~2| ) is inherited by h a p. Here 
and throughout the indices a,b, . . . run from . . . D and the indices a, f3, . . . will run 
from 2 . . . D and are the coordinates on B. 
A calculation then gives 

. . 1 , X ,\ 1 d2 u 1 d2 u 



+ 



-f g' 4 - d 



2/5 2# 2 gr 



d 4 

-i-hap nh aj 3, (14) 

dr gr z 
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where (A^/i)^ is the Lichnerowicz operator on B. All the other components of 
(Ai/i) a (, are zero because of the transverse tracefree property. This expression is the 
backbone of all the calculations in this paper. 

It should be noted that (13) are strong conditions in low dimensions. For the 
four dimensional Schwarzschild solution, for example, there are no perturbations of 



this form because S 2 does not admit any tensor harmonics [12]. We are looking at a 
potentially unstable mode that is specific to higher dimensional spacetimes. 

Our strategy in applying this to static spacetimes in the next section will be 
to calculate first a criterion for instability in terms of the minimum Lichnerowicz 
eigenvalue, X m in, on the base manifold B. In a later section we will then find this 
minimum for several relevant manifolds. 



2.2 Gauge Freedom 

Diffeomorphism invariance of the Einstein equations implies gauge invariance of the 
linear theory under 

Kb -> Kb + V a £ b + V b Ca- (15) 
The invariance can be used to set 

V a h ab = V a {Kb - \gabh c c ) = 0. (16) 

This is the transverse gauge condition which may always be imposed. The residual 
gauge freedom is given by vectors £ satisfying 

□6, + R a % = o. (17) 

Recall that the trace transforms as h a a — * h a a + 2V a £ a , we would like to find a £ 
satisfying (|i~7| ) and such that h a a — * 0. 

We show in appendix B that if the background spacetime is vacuum, possibly with 
a cosmological constant, then one may impose the transverse tracefree condition for 
perturbations as a gauge choice. This is slightly more subtle than the standard 
argument in which the cosmological constant is zero. 

Therefore the transverse tracefree choice made in the previous subsection is merely 
a gauge choice if the background spacetime is vaccum, with or without a cosmological 
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constant. However, this will not fix all the gauge freedom and we need to check that 
any solutions we find are not pure gauge. A pure gauge solution would be of the form 

Kb = V„6 + V fe £ a . (18) 

In appendix C we show that none of the modes considered in this paper are pure 
gauge. 

3 Application to static metrics 
3.1 Sturm-Liouville problem 

In this section we consider static metrics which solve the vacuum Einstein equations, 
possibly with a cosmological constant 

R ah = c{d + l)g ab . (19) 

This requires g = 1/f and the metric on B will be Einstein with 

R af3 = e(d - l)g af3 , (20) 

with e = il or e = 0. Tildes denote tensors on B. The cases e = ±1 correspond to 
having the same scalar curvature as S d or H d . The function / must be of the form 

f(r) = e-^Y~ 1 -cr 2 . (21) 

We look for unstable modes of the form 

h a ^x) = h af3 {x)r 2 ip(r)e ujt , (22) 

where x are coordinates on B and 

{A L h) a/3 = Xh aP . (23) 

As discussed in the previous section and in appendix A, we expect these tensor modes 
to be the dangerous modes. This is similar to the situation encountered in recent 
studies of stability of AdS p x M q metrics 




S 



The pertubation must satisfy 5R a p = c(d + l)h a p and this gives an equation for 
cp that may be cast in Sturm-Liouville form 

It is convenient to rewrite this as a Schrodinger equation by changing variables to 
Regge- Wheeler type coordinates and rescaling 

dr ± 
dr* = — , $ = rzip. (25) 



Equation (24) now becomes 

dr 2 

where the potential is 



+ F(r(r*))$ = = E$, (26) 



y ( r ) = 7^ + ^ — - + 4 72- 2c ( d + 1 )/- ( 27 ) 



Thus the stability problem reduces to the existence of bound states with E < of the 
Schrodinger equation with potential V(r). If such a bound state of the Schrodinger 



equation exists then the spacetime fllOD is unstable to modes of the form (22). That 
is to say, there will be an instability if the ground state eigenvalue, Eq, of (26) is 
negative. 

The normalisation of wavefunctions must take into account the weight function 



of (24), but the usual normalisation is recovered for $ 



1 = / ^ dr = l = / * 2 dr„. (28) 



f J f 

This condition of normalisability that is necessary to set up the Sturm-Liouville prob- 
lem is just the condition of finite energy of the gravitational perturbation (^). The 
background spacetimes ( [Tc| ) have a timelike Killing vector £° = 1, up to consideration 
of horizons. This allows the total energy of the perturbation on a spacelike hypersur- 
face with normal n° = 1/ f 1 ^ 2 to be well-defined, independently of the asymptotics 
of the background spacetime []24| 

E oc J t^n^g^)d d+l x = J jt 00 ^d d 3:dr, (29) 
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where too = G q = Rqq — ^gooR^- 2 \ the second order change in the Einstein tensor 
under the perturbation. Now note that the kinetic part of ( pS| ) contains terms like 
„d 



h af3 d d h a pyjg( d )d d xdr oc e 2ujt J h af5 h a ^ g( d )d d x J p 2 —dr (30) 



oc / ip -jdr. 



Thus requiring finite energy will recover the normalisation (28). 

Besides normalisation, we must also consider boundedness properties. The linear 
approximation to the equations of motion requires that h a b <C 1. For black hole 
spacetimes with event horizons, one should re-express solutions in Kruskal coordi- 



nates 25] near the horizon and check boundedness there. This is because Kruskal 
coordinates are well-behaved at the horizon. Fortunately, the mode we are con- 
sidering has no t or r components and therefore is essentially unchanged in Kruskal 
coordinates. Thus it is sufficient to check boundedness in r in the original coordinates 
of (10). However, in §3.2.2 we show explicitly boundedness in Kruskal coordinates 



for completeness. 

From (p2|), we see that boundedness requires p{r) = $(r)r~2 < 1. As r — ► oo 
the function $(r) must go as rs or lower power of r. This is a weaker constraint 
than is imposed by finite energy (|28|), with or without cosmological constant. As 
r^Owe must have <3?(r) going as r 2 or higher power of r. This will almost always 
be a stronger constraint than that required by finite energy (^). However, in many 
of the applications in this section there will be an event horizon at some finite ro, 
where /(ro) = 0. The condition of boundedness will then simply be that $(r) is 
bounded at ro. The finite energy condition will be that <£(r) goes to zero on the 
horizon, because the zero of f(r) is simple. In the cases below, we need to impose 
the stronger condition for each limit. However, in almost all cases we encounter, the 
soutions either satisfy both or neither of the criteria. 

3.2 Vanishing cosmological constant 
3.2.1 Asymptotic criterion for stability 

Set the cosmological constant c = 0. Asymptotically / — > 1 and r = r*. Alternatively, 
this is the massless case. We will derive first a criterion for instability by solving the 
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asymptotic Schrodinger equation ( pq ) with / = 1, and hence 

d 2 - lOd + 8 + 4A 
Mr) = — 2 , (31) 

and requiring suitable behaviour in the interior. Call the asymptotic solution ^oo* 
The range here is < r < oo. The argument of this subsection is in fact independent 
of the interior form of /. 

The asymptotic solution which decays at infinity is 

1 



* 0O (r) = K[riK v (ur)\, v = - y/(5 - d) 2 - 4(4 - A), (32) 



where K v {ojr) is the modified Bessel function that decays at infinity [26|. The be- 
haviour of ( |32"|) for small r and real, positive v is $oo( r ) ~ j--^ 1 / 2 . Three cases 
should be distinguished. If v > 1 the solution is divergent and not normalisable 
according to (|28|) . If 1 > v > 1/2, the solution is divergent but normalisable. If 
1/2 > v > 0, the solution goes to zero for small r. Another possibility is that the 
index v = iv^ is pure imaginary, in which case the Bessel function oscillates in the 
interior as sin(i/j lnr) and the wavefunction $ is then normalisable. We see that non- 
divergent normalisable solutions occur precisely when the potential (|3l| ) is negative, 
that is 

d 2 - lOd + 8 + 4A o 1 

1 =^ 2 -4^°- (33) 

We have finite energy solutions for a continuous range of lu > 0. The continuous 
spectrum of arbitrarily low energy is a direct consequence of the asymptotic potential 
(|3l|) being unbounded below. This will not be the case for the full potential and the 
spectrum will become discrete. 

In all of the cases of the previous paragraph, (p = $r~2 is not bound at the 
origin and so none of these solutions give instabilities of the massles, / = 1, metric. 



However, following [27| we note that in the oscillatory solutions with imaginary index, 
the derivative takes all values and so one might expect to be able to match the 
asymptotic solution to an interior solution for which / ^ 1, at least for certain 
discrete values of uj. Thus if a A exists such that v is imaginary then the Schrodinger 
equation should have a bound state and the metric is unstable. That is 

(5 — d) 2 

^min < A c = 4 - - — — — 44> instability. (34) 
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This is the criterion for instability of a massive black hole. We have also shown that 
the massless case is always stable. Concrete examples of Lichnerowicz spectra giving 
stable and unstable spacetimes are given in §5. 

The vacuum solution for / is of course (^Tj) which is now just the Schwarzschild- 
Tangherlini [13] black hole and the Asymptotically Conical (AC) variants considered 
in H |6| . The radial function is 



n 



d-X 



/W = l-(-J • (35) 

The higher dimensional Regge- Wheeler tortoise coordinate ( |25| ) may be given explic- 
itly in this case as [^] 

d-l 27rm/(d-l) . . 

r, = r + Y, — aln[r-e 2mn ^ d ^aj . (36) 

n=l 

There is an event horizon at r = a, and so the range of r is a < r. The potential 
becomes 

10d-8-4A /a^-id 2 ' 
r) 4 



1 /a\ d + 1 
V(r) = V OD (r) + -(^ 



(37) 



Oi 

It follows that V(a) = 0, as was clear from the initial definition (|27|), and that 
V(r) — * as r — > oo. 

The potential ( ^7|) can be seen to be always positive for a < r < oo if d 2 — Wd + 
8 + 4A > 0. This was the condition for the asymptotic potential to be positive also 
and signalled the nonexistence of finite energy solutions $ in this case. Thus, as 
expected, there are also no solutions to the full equations in this case. 

To establish the criterion for instability (|34|) we still need to check that there are 
solutions when the asymptotic solution oscillates in the interior and that there are 
none for the range 1/2 > v > 0, where an asymptotic solutions exists but does not 
oscillate. These statements will be supported numerically in the next section. The 
conclusion will be that a generalised black hole is unstable if and only if the base 
manifold has a Lichnerowicz spectrum satisfying $34). 

3.2.2 Numerical support for the asymptotic criterion 

The Schrodinger equation (26) with potential ([37]) has a regular singular point at 
the event horizon r = a. Thus we may perform a Taylor expansion of the solution 
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of the equation about this point. The leading order terms are found to be $ ~ 
(r — a) ±ctu '/( d ~ l \ so long as the exponent is non-integer. Typically there is thus 
one divergent and one convergent solution at the horizon. Further, the convergent 
solution vanishes on the horizon and therefore satisfies both the finite energy and 
boundedness requirements. We would like to see whether the solution that is well 
behaved at the horizon is also well behaved at infinity, giving a bound state. 

Before solving the equation we can check explicitly, following [28], that the well 



behaved solution remains well behaved at the horizon in Kruskal coordinates. In- 
cluding the time dependence and using the limit of (|36| ) as r — > a, we see that the 
mode (E2I) behaves near the horizon as 



Kp ~ (r - a^'^e^Kp = e" (i+a(r - Q)/( ^ 1)) /i a/3 ~ e" {t+T * ] h a(3 . (38) 
The Kruskal coordinates R, T are given by 

R±T ~ e /»(r*±*)/2 = e (d-l)(r„±t)/2a (39) 

Therefore the mode goes as 

h a p~(R + T)* u /f , Wh afl , (40) 

which is well behaved on the future horizon R — T = 0. It is not difficult to see 
that this expression will remain true for other functions f(r), such as that for the 
AdS black holes studied below. In this case / should of course be evaluated at the 
horizon, which would no longer be a. 

To investigate the equation numerically we first find a series expansion about 
the horizon of the solution that is regular at the horizon. We use this to set the 
initial conditions away from horizon itself. By taking sufficient terms in the series, 
this may be done to high accuracy. We then choose a (positive) value for uj and the 
equation can be numerically integrated. The solution will always diverge for large r 
because it is extremely unlikely that the ui we have specified corresponds precisely 
to a bound state. However, by varying uj we may see that the solution diverges to 
positive infinity for some values of u and to negative infinity for others. Because 
the solutions of the differential equation depend continuously on the parameters of 
the equation, there must be a bound state for some intermediate value of uj. The 
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solutions can be double checked by integrating in from infinity towards the horizon, 
although this is less accurate when cj is small, as for the interesting cases. Tables 1 
and 2 show values of uj between which the lowest lying negative energy bound state, 
u max = ~~ Eq, is found for various small values of d and A c — A, where A c is the critical 
value of A of (p4[). Without loss of generality we take a = 1. The dependence of the 
eigenvalues on a is determined on dimensional grounds to be u oc 1/a. Another way 
of seing this is that the Schrodinger equation is invariant under a — > ka, r — > kr, 
ui — * u)/k. 

Table 1: Lowest bound states for a = 1, d = 4 (D = 6) generalised Schwarzschild. 



A c -A 


Lower bound for uj max 


Upper bound for u max 


0.5 


3.63 x 10~ 2 


3.654 x 10" 2 


0.2 


2.64 x 10~ 3 


2.67 x 10~ 3 


0.1 


1.41 x 10" 4 


1.416 x 10" 4 


0.08 


4.341 x 10~ 5 


4.35 x 10~ 5 


0.05 


2.22 x 10~ 6 


2.28 x 10~ 6 


< 


No solutions found. 



Table 2: Lowest bound states for a = 1, d = 8 {D = 10) generalised Schwarzschild. 



A c -A 


Lower bound for ui max 


Upper bound for uj max 


0.5 


2.28 x 10~ 2 


2.29 x 10~ 2 


0.2 


1.59 x 10~ 3 


1.595 x 10~ 3 


0.1 


8.391 x 10~ 5 


8.399 x 10" 5 


0.08 


2.58 x 10~ 5 


2.59 x 10~ 5 


0.05 


1.34 x 10~ 6 


1.35 x 10~ 6 


< 


No solutions found. 



The results in Tables 1 and 2, and other similar results for different values of d, 
suggest that as A approaches the critical value from below, the energy of the lowest 
bound state rises and tends towards 0. This would provide a nice realisation of our 
expectation from the previous subsection that there should be no negative energy 
bound state if A > A c . Thus instability is according to the criterion (|34|). 
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3.3 Finite cosmological constant 
3.3.1 Topological black holes 

Let the cosmological constant c = —L 2 and let the base manifold B have negative 
curvature, e = — 1. The metric is seen to have no cosmological horizon. If a = it 
has no singularity and an event horizon at r = These are the so-called topological 
black holes 0, [l|, EI- Thus / is 

f[r) = -l + LV. (41) 

If B is a hyperbolic manifold B = H k /T, with T C SO(k, 1) a suitable discrete group, 
then the metric is locally Anti-de Sitter space. 
The potential (27) is 

~d 2 - 10rf + 8-4A L 2 (d 2 + 2d)' 



V(r) 



(1-LV). (42) 



4^2 4 

As expected, the potential vanishes on the horizon. The potential is not necessarily 
positive outside the event horizon, so a priori there exists the possibility of bound 
states with negative energy if — d 2 + Wd — 8 + 4A < 0. The Schrodinger equation ( p^ ) 
may be solved exactly in this case. Two solutions are 

$ ± ( r ) = r (l±C)/2 (1 _ L 2 r 2 ) -a,/2L x 

„ / u l-d±C uo 3 + d±C 2±C r2 2 \ lt s 
^{-2L + —^--2L + —i—— ;iV J' <" 3 > 

where C = yjjd - 5) 2 - 4(4 + A) and 2 Fi (a, 6; c; x) is the hypergeometric function. 
These are generically linearly independent. 

Because there is no cosmological horizon, the perturbation extends to r — > oo. To 
consider the asymptotics of ( |43l) , use the following result for hypergeometric functions 



2Fi(a,b;c;x) = kU-x) a 2 F\ [ a, a - c + 1; a - b + 1; - 

x 

+hi-xy b 2 Fi (6,6-c+l;6-a + l;i ) (II) 
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whereQ 

r(c)r(fe-a) 

1 T{c-a)T(by 
F(c)T(a-b) 

k 2 = — ■ (45) 

r(c — o)r(a) 

Using the series expansion of the hyper geometric function about the origin, this 
implies that as x — > oo 

2 Fi(«,6;c;x) = fci(-x)- a fl+ a(a " C ^ 1 1) - + 

\ a — o + 1 x 
- b ( , . Hb-c+1) 1 



In particular, the power law behaviour at infinity is x max (- 5R ( a )>- !R ( 6 )) . j n both the 
solutions of (p3|) we have — 3?(a) > — 3? (6), so the x~ a term is dominant as x — > oo. It 
is then easy to check that the overall leading asymptotic term of both the solutions 
<3?± at infinity is C(r rf//2 ). These solutions are never normalisable in the sense of (28). 



However, because infinity is a regular singular point of this Schrodinger equation, it 
is possible to take a linear combination of these two solutions that gives the other 
allowed power law asymptotics, in this case 0(r~ ( 2+d )/ 2 ). This is precisely the linear 
combination of <f>± in which the x~ a terms of ([44] ) cancel. The result is, using the 
symmetry 2-^1(0, b; c; x) = 2Fi(b,a;c;x), 

$ 3 ( r ) = r -(a+«0/a+«/^(i _ lV 2 )-^ x 

/ to 3 + d + C uo 3 + d-C6 + 2d 1 \ 



2L 4 ' 2L 4 ' 4 ' L 2 r 2 , 

This solution has acceptable behaviour as r -» 00. We now need to check the be- 
haviour as r —* -f. To do this we use another identity of hyper geometric functions 



2-Fi(o, b; c; x) = /ii 2-^1(0, a + & + 1 — c; 1 — 2) 

+ /i 2 (l - a;)^ a ^ 2 Fi(c - a, c - 6; 1 + c - a - 6; 1 - x), (48) 



1 There is a subtlety here which is that if d is odd then a — b is a negative integer and one of these 
gamma functions diverges. However, the solution (^7]) that we obtain is the solution which decays 
at infinity even in these cases. 
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where 

T(c)T(c-a-b) 



hi 



T(c-a)T(c-by 



2 T(a)T(b) ■ 1 j 

Applying this to (^7j) we see that generically for the modes with uj > that we 
are looking for, the leading term as r — ► \ is O ((1 — L 2 r 2 ) _aj / 2Z/ ) . Inserting this 
and ( pi] ) into the normalisation condition (|2^) we see that the solution is not square 
integrable at the horizon. Furthermore, the solution does not satisfy the boundedness 
requirement because $, and hence <p, diver gs at the horizon. However, if C is real 
and C > 3 + d for some mode, which requires 

Amm < -4d, (50) 

then we may set 

- = (^f^) L > 0. (51) 



For this mode we have that hi = and h% = 1 in ( |48| ) because one of the T functions in 
the denominator of h\ diverges. Another way of seing this is that the hypergeometric 
function is just a polynomial in this case. It then follows from (|48| ) that the mode now 
has a better behaviour at the horizon, going as O (l — L 2 r 2 ) u} / 2L ) . In particular it is 
normalisable and bounded at the horizon. Therefore it gives an instability. Because 
the base manifold has negative curvature, the arguments of Appendix A do not apply 
in this case and therefore we don't know about the effect of scalar and vector modes. 
The statement we can make is that if the base of a massless topological black hole 
has a Lichnerowicz spectrum satisfying $5(\ ) then it is unstable. Some results on the 
Lichnerowicz spectrum for negative scalar curvature Einstein manifolds are collected 
in §5. 

3.3.2 Brane world metric 

Here f(r) = r 2 . The base manifold is Ricci flat because e = 0. The metric is more 
familiar in terms of z 



r 



ds 2 = [alz 2 - dt 2 + dsfj . (52) 
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The r coordinate extends from the origin to infinity. This metric is of the general 
brane world form considered in (I?], [l^]. More general such metrics will be studied 
in the section below on time dependent solutions. 
The potential ( p7| ) is 

V(r) = X + ^+f> r \ (53) 

The starred coordinate (25) is just r* = —z. The Schrodinger equation ( p6| ) may be 
solved exactly in this case. The general solution is 

*(r) = Ar^% +d)/2 (^±^) + Br-^K (1+d)/2 (^±^), (54) 

where A, B are constants and I U ,K U are the modified Bessel functions. We need to 
take the real or imaginary part depending on whether VA + uj 2 is real or imaginary. 
The term with the l v function is always normalisable (pq) as r — ► oo, going as 
Q(j.—(2+d)/2^ ^he ^ erm w [t\! the K u function diverges as as r — > oo and 

hence is never normalisable because d > 1, although ip is bounded. We must then 
check the I v solution as r — ► 0. In terms of z, the normalisability condition, with an 
implied real or imaginary part being taken, is 

/•oo 

J z lfi+d)/2( z V / ^ + oj 2 )dz < oo. (55) 

If A + to 2 > then the integrand diverges exponentially as z — > oo. There was no 
chance of a solution in this case because the energy, —to 2 , would have been lower 
than the minimum of the potential A. If X + uj 2 < then the integrand oscillates with 
constant magnitude as z — > oo. In either case, normalisability cannot be enforced. 
Furthermore, (p = <£r~2 is seen to diverge as r — ► and so the boundedness condition 
is not satisfied. Therefore there are no finite energy solutions and the brane world 
metric (f5lL) is stable under this perturbation, independent of the base manifold. 



3.3.3 Schwarzschild-Tangherlini-Anti-de Sitter black holes 

Set c = —L 2 and e = +l. Ifa = then we have a generalised Anti-de Sitter space 
that describes the asymptotics of a Schwarzschild-Tangherlini-Anti-de Sitter black 
hole. The function /(r) is 

f{r) = 1 + LV. (56) 
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There are no horizons. The potential is now 



^oo(r) 



d 2 - lOd + 8 + 4A {d 2 + 2d)L 



(L 2 r 2 + 1). (57) 



4r 2 4 
Two exact solutions to the Schrodinger equation are 

* ± ( r ) = r (l±C')/2( 1 + L 2 r 2 ) -to/2I x 

„ f-iui l-d±C -iuj 3 + d±C 2±C 

2-F1 1 , 1 ; ; -£V , (58) 



with C" = yj{d — 5) 2 — 4(4 — A). This is fairly similar to the topological black hole 
case, but we now need to consider different limits. By exactly the same arguments 
as for the topological black hole, the solution that is well behaved at infinity, going 

as 0(r-( 2+d )/ 2 ), is 

$ 3 ( r ) = r -C2+d)/2+iu J /L( 1 + L 2 r 2ytu J /2L x 

„ f-ioo 3 + d+C -iu 3 + d-C 6 + 2d -1 \ . s 
2i?1 (^ + ^'^ + ^^wJ' (59) 

By using (|44]) again, we see that the behaviour as r — > is O ^r^^ c )/ 2 ^ . This both 
converges and is normalisable if d 2 — Wd + 8 + 4A m j n < which is unsurprisingly 
also the condition for the potential (57) to be negative and indeed unbounded below. 
There is a continuum of negative energy bound states with u > in this case. 

_ d 

None of these solutions satisfy the boundedness condition because cp = <5r 2 always 
diverges at the origin. Thus the massless case is stable against the perturbation. The 
phenomenon of the hypergeoemtric series terminating for special values of u to give 
a well behaved mode does not occur here because of the i in front of the u> in 
More interestingly, if 

A min <4-^=^, (60) 

then the solution is oscillatory in the inner regions. The oscillatory behaviour suggests 
we can make a statement about massive black holes also. The condition is clearly the 
same ([34|) as we found before in the case of vanishing cosmological constant. This 
is perhaps not surprising given that the potentials, (|3l|) and (|57|), are the same near 
the origin. The comments of §3.2 should go through. However, there is now a second 
lengthscale in the problem, 1/L, which could delay the onset of the oscillations until 
beyond the event horizon, in which case there will be no solution and no instability. 
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A few things may be said more concretely. In the massive case 

'a 



./ID (~J ; • /-V. (61) 



The potential becomes 



V(r) = Voo (r) + - 2 (-) (-) 



lOd - 8 - 4A - 2dLV /a\«^ 

T 



(62) 



The potential has the expected property that it vanishes at the horizon, where fir) = 
0. Furthermore, when d 2 — lOd + 8 + 4A > it is everywhere positive, as was the 



asymptotic (massless) potential (57). Thus there will be no instability in these cases. 
We now need to see numerically when a solution exists and what the role of the new 
lengthscale is. 

3.3.4 Numerical results for Schwarzschild-Tangherlini-Anti-de Sitter 

We wish to use the methods of §3.2.2 to examine the effect of the new lengthscale 
1/L. First note that the Schrodinger equation is now invariant under the scaling: 
a — > ka,r — » kr, uj —* oj/k,L — > L/k. Previously, when there was no L we used this 
to set a = 1, which was the location of the horizon. Again, we want to scale the 
horizon to 1. This will now require scaling so that a = (1 + L 2 ) 1 /^ -1 ). There is a 
scale-invariant dimensionless mass 

M = aL, (63) 

which allows us to talk about large and small black holes independently of the scaling 
used. We expect the criterion for instability to be the same as for Schwarzschild black 



holes (34) when the AdS black hole is small M — > 0. As we increase M we expect 
the black hole to be stabilised by the cosmological constant. This behaviour indeed 
happens and is illustrated in table 3. The numerics were done as in §3.2.2 For the 
case d = 4 (hence spacetime dimension D = 6) and with A satisfying the instability 
criterion by A = A c — 1 we see that the unstable mode is stabilised if M > 0.058.... 

This will be the generic behaviour. A given base space will have a minimum 
Lichnerowicz eigenvalue, A. // this is less than the critical value A c then we can find 
a critical value for the dimensionless mass, M c , such that if M c < M then the unstable 
mode is stabilised. If M c > M then the AdS black hole is unstable. Alternatively, we 
could think of the mass as altering the expression for A c . 
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Table 3: Lowest bound states for A = A c — l,d = 4 (D = 6) generalised 



Schwarzsc 


tiild-AdS. 


L 


1V1 = CtL 


Lower bound for u> max 


Upper bound for ui max 








0.1331 


0.1335 


0.01 


0.0100003... 


0.131 


0.132 


0.05 


0.04004... 


0.0591 


0.00599 


0.055 


0.05505... 


0.0264 


0.0269 


0.057 


0.05706... 


0.0005 


0.00051 


> 0.058 


0.05806... 


No solutions found. 



3.3.5 Schwarzschild-Tangherlini-de Sitter black holes 

Set c = L 2 and e = +1. This gives us a generalised Schwarzschild-Tangherlini-de 
Sitter black hole. We have, 

/(r) = l-(^)"-LV. (64) 

There is a cosmological horizon at finite radius, and therefore we cannot discuss 
an asymptotic solution to the Schrodinger equation, because the mass term is not 
negligible near the horizon. We have all the information necessary to tackle this 
problem numerically, but various cases must be considered separately depending on 
the values of a and L. This is somewhat out of the main line of development of this 
work and so it will not be considered here. 



4 Time-dependent solutions 

4.1 Generalised de Sitter 

The metric form (10) also covers a range of cosmological solutions. For example, a 
generalised de Sitter metric may be written 

ds2 = L 2 r t_ ! + r2d ~ s ^ ( 65 ) 

where r is the time coordinate now and L is a constant. This is of the form (|1C 
with / = 0. Note that d would now be 3, not 2, in the usual four dimensional case. 
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Consider a perturbation, using the same notation as in (f22[), 

h a/3 (x) = h a/ 3(x)r 2 (p(r), (66) 

and impose the Einstein equation, 5R a p = dL 2 h a p. There is now a cosmological term 
which is that of ordinary de Sitter space in D = d + 1 dimensions. The equation for 
the perturbation <p is most familiar in terms of the coordinates 

cosh Li = Lr, (67) 

where the metric becomes 

ds 2 = _ dt 2 + ^J± d ^ (68) 

and the equation (|l~4| ) becomes the equation for a scalar field on de Sitter space 

d 2 ip dip L 2 v , , 

— + dLtanhLt-^- + « — (A + 2 - 2d )u> = 0. 69) 

dt 2 dt cosh 2 Lt 

This provides a check on our expression (14) and also shows that at late times as 
t — > oo the leading term in each of the two linearly independent solutions is 

<p ~ A + Be~ dLt . (70) 

With A, B constants. Thus perturbations are frozen in, independent of the dimension 
and the form of the base Einstein manifold B. This is just the behaviour of such 
perturbations in standard four dimensional inflationary metrics [pSfl. 

4.2 Generalised Anti-de Sitter 

Generalised Anti-de Sitter space can be treated similarly. Write the metric as 

^ 2 = n$Vr+^ ( 71 ) 

where again r is the time coordinate and L a constant. The base manifold B must 
now have negative curvature and would be H d for Anti-de Sitter itself. Consider the 
perturbation 

h a p(x) = h a/3 (x)r 2 ip(r), (72) 
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and impose the Einstein equation, 5R a p = —dL 2 h a p. This is as for the de Sitter 
case considered previously but with a negative cosmological constant. The familiar 
coordinates for the space are 

sin Lt = Lr. (73) 

The metric in these coordinates is 

ds * = _ dt 2 + s j^Jd dS l (74) 

These coordinates make explicit the Big Bang and Big Crunch singularities at t = 
and t = J, unless the spacetime is Anti-de Sitter. The equation for the perturbation 
becomes 

4^ + dLcotLt^f- + — % — (A + 2d - 2)<p = 0. (75) 
dt 2 dt sin 2 Lt 

The general solution to this equation is 

<p = A (sin Lt) (1 - d)/2 P (d _i)/ 2 C/2 (cos Lt) 

+ B (sin Ltf- d)/2 Q {d _ 1)/2 c '/ 2 (cos Lt), (76) 



where A, B are constants and as before C = y{d — 5) 2 — 4(4 + A). P u ^ and Q u ^ are 
Legendre functions of the first and second kind, respectively. These may be expressed 
in terms of hyper geometric functions as follows J30| 



Q„*(x) oc (^-i r /^-S Fl (^^|±i,f^±l ; , + ? ; J_). (77) 

In the present case a; = cos Lt with < t < J. Thus we need to check regularity 
properties at t = 0, ?, corresponding to x = 1, 0, —1. 

It is easy to check using ( |44| ) that both solutions are finite at t = that is x = 0. 
Further, it is clear from ([77]) and (|7^) that behaviour as t — > ^ will be the same as 
for t — > 0, up to phases in front of each Legendre function. In particular this means 



that regularity properties will be the same at these points. Both solutions of (76) 
diverge as C?(^ 1_c(_c '^ //2 ) as t —> 0. Therefore they also diverge as t — > ^. However, 



i = is a regular singular point of the equation (|75|) and therefore there will be a 
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linear combination of these solutions that has the other allowed power law behaviour 
as t — > 0, namely 0{t^ l ^ d+c ^ 2 ). This will be oscillatory and divergent if C is pure 
imaginary. If C is real it will converge if C > d — 1 and diverge otherwise. There is 
also the possibility that well behaved modes will exist for special values of C where 
the hypergeometric function becomes a polynomial. 

The main conculsion of the previous paragraph is that there are always modes 
that if excited at some finite time will diverge in the future. Thus this AdS cosmology 
is unstable, independent of the base manifold and the dimension. 

4.3 Ricci flat Lorentzian base (Brane world metrics II) 

Let the base 6 be a (i-dimensional Ricci flat spacetime. Suppose we know the spec- 
trum of Lichnerowicz modes on B with eigenvalues A and such that the modes grow 
in time. In particular, if this spectrum includes a zero mode, then the spacetime B 
is unstable. 

The spacetime B may be embedded in a D = d+ 1 dimensional Einstein manifold 
with negative scalar curvature 



ds 2 = \ (dz 2 + dsj) =%- + r 2 dsj, (78) 



For example, if B is a black hole metric we obtain a black string in AdS spacetime. 
The change of variables in (f7£|) is of course r = 1/z. We would like to see whether 
this spacetime is unstable under any of the growing modes in B. Consider the 
perturbation 

Kp(x) = r 2 ip(r)h a/3 (x), (79) 

where h a p is a Lichnerowicz eigenmode on the base with eigenvalue A. From (|l~4]), 
ignoring the terms with /s, the equation for (p(r) coming from A^hag = —2dh a p is 

^ + d + i^_V =0 _ (80) 

ar z r dr 
The general solution to this equation for A > is 

p(r) = Ar~ d / 2 I d/2 + Br- d l 2 K d/2 ( , (81) 
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where A, B are constants and Iu, Kn are the modified Bessel functions. If A < then 
the expression is most transparent if we let A — > —A and replace Iw, by the Bessel 
functions J^,Y^. Finally, if A = then the solution is 

tp(r) = A + Br~ d . (82) 

To see which, if any, of these solutions are acceptable we need to find the energy 
of the perturbations. Assuming that B has a timelike Killing vector, this is similar 
to the argument in §3.1 

E oc J t^n^ v ^g~^>d d x = j t^n^sJgV-Vr^drd^x 

~ J ip 2 r d - 3 dr, (83) 

where we used the metric (fT§|) and no oc r, £o oc r 2 , too ~ 9 J ( r ) 2 - Note that g^ 1 ^ is 
the spatial metric on B. It is easy to see that the normalisation condition (|83|) is the 
same as the norm of the equation (^) cast in Sturm-Liouville form. The eigenvalue 
of the Sturm-Liouville problem is now A and the weight function is seen to be r d ~ 3 , 
consistent with (^). The boundedness condition is this context is simply that <p{r) 
remains finite in the range < r < oo. 

If there is a A = mode, the d dimensional spacetime B is unstable. However, 
from (|82| ) it is clear that none of these solutions are normalisable in the d + 1 di- 
mensional spacetime. Further, all except the constant solutions are not bounded. 
Therefore the unstable mode does not carry over to the full spacetime because it no 
longer has finite energy. 

For A > we see that the term with the 1^ Bessel function has an exponential 
divergence as r — > whilst the term with the Bessel function goes to zero expo- 
nentially and hence is normalisable at the origin. However, this term goes as O(r ) as 
r — ► oo and hence, although bounded, is not normalisable at infinity because d > 2. 
Thus there are no normalisable solutions with A > 0. 

For A < 0, as r — ► oo, the term with the Bessel function goes as 0(r ), 
and therefore is not normalisable at inifinity. The term with the Bessel function 
goes as 0(r~ d ) and therefore is normalisable. As r — * 0, the term oscillates as 
0(A 1 ~ d ^ 2 cos - ) and hence are not normalisable ( |83| ) or bounded as r — ► 0. 
Thus there are no normalisable modes of this type. 
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In conclusion, embedding a Ricci flat spacetime B into a higher dimensional space- 
time with cosmological constant as in ( |78|) will stabilise at the linear level any unstable 
modes of B. Furthermore, no other unstable modes of the form we consider appear. 
This should be contrasted with a similar embedding into a higher dimensional Ricci 
flat spacetime where the stability properties get worse due to negative Lichnerow- 



icz modes in the initial spacetime [31], such as the Gregory-Laflamme instability of 



nonextremal black strings p2"| , |33| , 25]. It was argued in 34] that a Gregory-Laflamme 
instability existed also for black strings in AdS spacetime. However, the perturbed 
mode presented there, which agrees as a special case with the modes we have just 
considered, did not have h a p bounded as is required in a linearised analysis. The 
phenomenon of perturbations to brane world metrics diverging in the bulk has been 
observed before [^] and is related to the bad behaviour of the curvature at the hori- 
zon. Bounded modes, and hence the instability, could reappear if one modifies the 
setup, such as by adding a negative tension brane at finite position. 

One might worry that the instability of a very thin short black string should 
not be affected by immersion in an Anti-de Sitter spacetime with large radius. This 
may well be true. However, if one assumes that the black string runs all the way to 
the horizon there seems to be no way of avoiding our conclusions, although there is 
already a singularity near the horizon in the unperturbed metric p6|]. The methods 
used here can say very little about what would happen for a "cigar-like" black string 
configuration. 

4.4 Double analytic continuation 

Higher dimensional versions of the Schwarzschild bubble solution IS] have been con- 



sidered recently in a search for well-behaved time dependent backgrounds in which 



to study string theory [37]. Ultimately the Kerr solutions turn out to be more ap- 
propriate, but the Schwarzschild case contains many of the relevant features. The 
solution, obtained via double analytic continuation of the Schwarzschild- Tangherlini 
solution is 



r 



dr 2 

d^ 2 + -j— r + r 2 {-dr 2 + cosh 2 rdOL,), (84) 

1- {f) 
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where r > a, tp has period j^t and dftj_ l is the round metric on S l . The coor- 
dinate tp = it is the Wick rotated time from the black hole solution and it = 6 — ^, 
where # was the usual angular parameter on S d and r is now the time on dS^. Re- 
cently, such "bubbles of nothing" have also been considered by analytically continuing 



AdS-Schwarzschild black holes ||, ||. 

A classically stable background is needed for string theory. Classical instabilities 
are manifested in the string worldsheet theory as a renormalisation group flow arising 
from higher string loops with a fixed point that is not close to the original background 
in string units. 

It was argued in [37] that classical stability of (84) follows from the classical 
stability of the corresponding black hole solution. Any mode on the bubble must 
be periodic in the ip variable. But e~ tLU ^ = e wt , so these correspond to growing 
modes of the black hole with certain frequencies lo allowed by the periodicity of tp. 
Furthermore, the part of the black hole mode that is a harmonic on S d becomes a 
harmonic on dS^ with the same eigenvalue. In the four dimensional case of ||, ||, |l(]] 
these modes are scalar and vector harmonics, whilst in our case we are considering 



tensor harmonics 12]. The equation for the radial dependence is the same in the 
black hole and bubble cases and therefore the criterion for the existence of solutions 
is the same. On the sphere, the harmonics are trigonometric functions of 8 and these 
become hyperbolic functions of r, which correspond to growing and hence unstable 
modes on the bubble. 

The conclusion of the previous paragraph is that any unstable mode on the black 
hole with appropriate u signals an unstable mode on the bubble. Conversely, any 
mode on the bubble signals an unstable mode on the black hole. Thus we have 
extended the stability arguments of the higher dimensional Schwarzschild bubble in 
|37j] to include the tensor mode that is being considered throughout this paper. 

There is a subtlety, however. After doing the double analytic continuation, we 
must redo the calculations of energy because the timelike Killing vector has changed, 
if indeed there is still a timelike Killing vector at all. The calculation is very similar to 
the Lorentzian base case (|83l). There will typically be horizons on the base manifold 
across which the Killing vector changes sign, as in the de Sitter base of (|84|). In this 
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case, the integration over a spacelike hypersurface of the base should be restricted to 
the region inside the horizon. The integration over r is not changed. 

~ J ip 2 r d - 2 dr = J $ 2 r- 2 dr, (85) 

which is indeed different to the energy of the corresponding black hole perturbation 
(p9|). As a check of this expression, we can see that this is the same normalisability 
condition that we get from the Sturm-Liouville problem for ip. The Sturm-Liouville 
equation remains (24), but the eigenvalue associated with time evolution is no longer 



uj, but rather A, as was commented for the case of the Lorentzian base. The weight 
function is thus seen to be r d ~ 2 rather than the r d / f that we had for the black hole. 
This is in agreement with (p5|). The condition for boundedness remains the same, <p 
must not diverge. We may now go through the results from the black hole spacetimes 
of §3 again and recheck for finite energy using (|85|). We see that the existence and 
nonexistence of finite energy solutions remains the same in each case. 

If the Euclidean metric on the base manifold B of the generalised black hole 
solutions, (|lOl) and (]2l|), admits an analytic continuation to a Lorentzian metric, 
then one may construct a generalised "bubble of nothing" with metric 

dr 2 

ds 2 = f{r)d^ 2 + ^- )+ r 2 dsl (86) 

with the notation of (|84] ) and now ds^ is a Lorentzian metric. 

Assuming that the tensor harmonics on the Lorentzian base manifold are growing 
modes in time, the argument for the S d above case goes through also in the generalised 
case. Stability is thus related to the stability of a generalised black hole solution. We 
could also see this directly from the relationship between the Lichnerowicz operator 
on the full manifold and on the base (|l4|): perturbations of ( |86| ) that are periodic 
in ifi will give the same radial equation as perturbations of the generalised black 
hole metrics with exponential growth in t. Thus, for example, in the generalised 
Schwarzschild-Tangherlini case, the criterion (|34| ) will be the same. However, there 
is an important caveat: the frequencies that cause instabilities for the generalised 
black hole will only cause an instability in the bubble spacetime if they respect the 
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periodicity of ■0, that is 



u= N «- 1 \ NeZ. (87) 
2a 



In practice, this is unlikely to be true for any of the bound states in the discrete 
spectrum. Therefore we expect that the bubble spacetimes will generally be stable 
against the perturbation even when they violate the criterion (|34]). 

5 Spectra of the Lichnerowicz operator 

This section collects and applies mostly known results about the spectrum of the Lich- 
nerowicz operator on various Einstein manifolds. In our considerations of the static 
metrics above we found that in the case of Schwarzschild-Tangherlini, Schwarzschild- 
Tangherlini-Anti-de Sitter and topological black holes, stability of the spacetime 
against the perturbation (|i~2|), (13) depends on the spectrum of the Lichernowicz 



operator on the base manifold B (p4|). The objective of this section is to check that 
the cases we expect to be stable, i.e. round spheres, are stable and also to find exam- 
ples of spaces that are not stable, thus showing that the criterion for instability is not 
vacuous. It would be nice to have more spectra at hand, in particular for Einstein 
manifolds that are topological spheres such as the Bohm metrics H. We discuss 
first Einstein manifolds of positive scalar curvature, relevant for the Schwarzschild 
and Schwarzchild-AdS cases. We then discuss manifolds of negative scalar curvature, 
which are relevant for the topological black holes. 

In this section, tildes will be omitted from the metric perturbations on the base 
manifold B and we will consider transverse tracefree perturbations 

9<xf3 -> ga/3 + h a/ 3, h a a = V a /l a/3 = 0. (88) 

Define A by R af3 = Ag af3 . 
5.1 Round spheres 

When the base manifold is a unit sphere, Asymptotically Conical (AC) just means 
Asymptotically Flat (AF). The unit sphere B = S d has curvature 

Ra/3f8 = 9a-y9/3S ~ 9a89/3^: (89) 
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and therefore the Lichnerowicz operator is 

(A L h)ap = (-V 2 /i) Q/3 + 2dh aP . (90) 

However, the Laplacian on the sphere is a positive operator and therefore A m ; n > 2d. 
In fact, the spectrum of transverse tracefree tensor harmonics on S d , d > 3, is known 
to be @ 

(-V 2 h) af3 = [k(k + d-l)-2]h a/ 3. (91) 



But 2d is larger than the critical eigenvalue of (|34|) for all d, and therefore the 
spacetime is stable against this perturbation. In particular this means that the AF 
Schwarzschild-Tangherlini black holes are stable to the perturbation, as we should 
expect. The spectrum of quotients of the sphere by a finite group has the same lower 
bound and so the resulting ALF spaces also give stable black holes. 

5.2 Product metrics 

Let the base be metrically a product of Einstein manifolds. It is well known [41, ^IJ] 
that there is a Lichnerowicz transverse tracefree zero mode in which one component 
expands and the other shrinks, keeping the total volume constant. Specifically, if the 
metric decomposes as 

ds 2 d = ds 2 n + ds 2 d _ n , (92) 

then the mode 

Kp = ( J , (93) 

is easily seen to be tracefree, transverse and with zero Lichnerowicz eigenvalue. Com- 
paring with (|34|), this implies that the spactime is unstable with a product base for 
d < 9. It is intriguing that this is the same critical dimension for instability as was 
found in the context of generalised Freund- Rubin compactifications to AdS x B d pi]] , 
where the general dependence on the dimension is different. 
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5.3 Five and seven dimensional base 



In five dimensions, one has the Einstein manifolds T pq which are £7(1) bundles over 
S 2 x S 2 . See e.g. for the metric and other details. We will not make their choice 
of normalisation a = 1 and a will appear as a parameter related to the overall scale 
of the metric. Instead we require that A = 4, as this is the normalisation of the unit 



sphere in five dimensions. We have [22 

2 



(4a 2 - l)(12a 2 - l) 2 ' 



which is seen to imply 



1 



< a z < 



The minimum Lichnerowicz eigenvalue is 

Amir?. 



12a 2 - \/784a 4 - 240a 2 + 20 



1 



(94) 



(95) 



(96) 



The critical value in four dimensions (34) is A c = 4. For the range of a in (p5|), we see 
that \ m in < A c , with equality ocurring for the T 11 case where a 2 = Thus the T 11 
black hole is stable whilst all the other T pq black holes are unstable. This is precisely 
the behaviour that was observed in the context of Freund-Rubin compactifications 
in (22]. 



In seven dimensions, the Einstein manifolds M pqr , which are U(l) bundles over 
l2 x S 2 , and Q nin2res , which are U(l) bundles over S 2 x S 2 x S* 2 , have been studied 



in the context of Freund-Rubin compactifications [42, 43 1. Again, results will be 
quoted. Interestingly enough, in seven dimensions the manifolds with the required 
normalisation, A = 6, turn out to be stable precisely when the corresponding Freund- 
Rubin supergravity compactifications are stable. Thus the bounds we obtain are 
familiar. 

For M pqr , the stability depends on - and A as follows 



X, 



where x is defined by 



A 



4 + 4x- 2(25 -48x + 32x 2 ); 
l + 2x 



2x - 1 
x 2 (3 - 2s)" 



(97) 



(98) 
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This implies \ < x < |. We are interested in the case A = 6. Given that the critical 
value in seven dimensions is A c = 3 (|3~4|), we have that the solution will be stable to 
this perturbation for 

9 39 

n *0.64...<s<1.15...«-, (99) 
and unstable for x outside this interval. Translated into an interval for - this becomes 



7y/E 

0.63... < 



P 



17V66 . , 

<1.1 8 ...~_X_. (ioo) 
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Thus for this range of p, q we have solutions that are stable against perturbations of 
the kind considered in this work, and therefore stable overall in the cases of appendix 
A. This means that the Schwarzchild-Tangherlini solution with the sphere as base 
space is not the only classically stable solution of the form (|i~0|), contrary perhaps 
to initial expectations. These spaces are not necessarily simply connected, having 
fundamental group Z r . 

The situation is more complicated for Q n ^ n ^ n 3 . "We have \ m in = 7mmA, where 
'jmin is the smallest root of 

7 3 — 67 2 + 20(aiQ2 + CK2«3 + «3Cki)7 — 56aia2«3 = 0, (101) 
with the ajS defined by 

ni\ 2 ai(l + a 2 ) 2 



i - /-I i x 2 + c y lic 

n 2 J a 2 (l + ai) 
a 1 + a 2 + a 3 = l. (102) 

We will look at the simplified case in which a 2 = «3, implying n 2 = n.3. Further, set 



A = 6. The resulting range of stability to this perturbation is found to be [43] 



ss 0.85... < 



ri2 

m 



17V66 . , 

<1.1 8 ...~_X_ (103) 



Again we have found a countable infinity of AC Schwarzschild-Tangherlini solutions 
that are stable to this perturbation. The upper bound is the same as previously. The 
fundamental group is Zfc, where k is the greatest common divisor of 711,712,713. 

Finally, because the bound for stability on the Lichnerowicz operator is the same 
for a seven dimensional base as the Freund-Rubin bound, all the manifolds that 
are super symmetric in the supergravity context give a stable spacetime. A list 
of such spaces can be found in pfl ; examples are the squashed seven sphere and 
SO(5)/SO(3) max . 
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5.4 Negative scalar curvature manifolds 

Hyperbolic space and its quotients by appropriate finite groups have curvature 



RafS-yS = -9^9/36 + 9a59p~f> (104) 

and therefore the Lichnerowicz operator is 

(A L h) aP = (-V 2 h) af3 - 2dh a p. (105) 

The rough Laplacian is a positive definite operator on normalisable modes and there- 
fore the Lichnerowicz spectrum is bounded below by —2d. Comparing with (|50|) we 
see that the corresponding topological black holes are always stable against the tensor 
perturbations. 

The spectra of negative scalar curvature Einstein-Kahler manifolds was studied in 
p5| with the result that for these manifolds the Lichnerowicz operator on transverse 
tracefree tensors is bounded below by 2 — 2d. Therefore these manifolds always give 
black holes that are stable against such perturbations. 

6 Conclusions and open issues 

Let us summarise the results of this paper. 

(a) Spacetimes of the form ([n]) admit perturbations that are transverse tracefree 
tensors on the base Einstein manifold. The stability of these spacetimes against such 
perturbations depends on the spectrum of the Lichnerowicz operator on the base 
manifold. 

(b) Higher dimensional generalised Schwarzschild black holes are stable if and only 
if the Lichnerowicz spectrum on the base manifold is bounded below by the critical 



value of (34). This statement includes all perturbations to the spacetime. The same 
criterion holds for Schwarzschild-AdS black holes, although in this case large black 
holes become stabilised. Examples of base manifolds leading to stable and unstable 
black holes are given in §5. 
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(c) (Generalised) Topological black holes have a criterion for instability due to tensor 



modes on the base (50). An Einstein-Kahler base always gives a stable black hole, as 
do quotients of hyperbolic space. The brane world metrics of §3.3.2 are always stable 
against these tensor perturbations. Time-dependent tensor pertubations in gener- 
alised de Sitter spacetime, considered as a cosmological spacetime, are always frozen 
in whilst in generalised Anti-de Sitter spacetime they always include an unstable 
mode. 

(d) A spacetime with no cosmological constant may be embedded in a higher dimen- 
sional spacetime with a negative cosmological constant. Lichnerowicz zero modes 
which correspond to instabilities of the initial spacetime become stabilised. At the 
linearised level, no new instabilities of the form we consider are introduced. In partic- 
ular this means that the perturbative Gregory-Laflamme instability does not occur for 
AdS black strings in five dimensions, contrary to previous claims [34]. See discussion 
in §4.3. 

(e) Double analytic continuation of generalised black hole spacetimes, if admiss- 
able, produces a generalised "bubble of nothing" spacetime. These spacetimes are 
generically more stable to the analytically continued tensor perturbation than the 
corresponding black hole. However, one needs to redo energy calculations because 
the time direction has changed. 

Some issues remain to be addressed. For the perturbation analysis, these include 
a numerical study of the stability of generalised de Sitter black holes, an explicit 
calculation of the vector and scalar modes to confirm the arguments of appendix A 
and to extend them to cases where the base manifold has negative or zero curvature 
and an extension of the arguments to non-vacuum solutions. 

Regarding the Lichnerowicz spectrum, the spectrum is not known for many in- 
teresting metrics, such as the Bohm metrics. It is possible that the methods of the 
present work may be used to study aspects of this spectrum. 

In other directions, it might be interesting to see whether rotating black holes 
admit a similar generalisation. Also, to see whether these instabilities are useful in 
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the context of gravity- gauge theory dualities. 
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A Scalar and vector modes 

We argue here that if the Einstein base manifold B is compact, Riemannian and with 
e = 1, then the scalar and vector modes will not produce instabilities in the space- 
times under consideration. In particular, this covers the generalised Schwarzschild- 
Tangherlini(-AdS) cases in which we found a criterion for instability from the tensor 
modes. 



The argument has two steps. Firstly, we will recall [21] that the scalar and 
vector second order differential operators on B are bounded below by their minimum 
eigenvalue on the sphere S . Secondly, we will show that the equations of motion for 
these modes depend on the base manifold only through these differential operators. 
We will use h and h a to denote generic scalar and vector modes on the base manifold. 

The second order differential operator on scalars is Agh = — V Q V a /i. This is 
non-negative on compact manifolds without boundary. There is always a zero mode 
corresponding to a constant field. In particular, zero is the minimum eigenvalue on 
S d and also on any other such manifold B. 

The second order differential operator on vectors in B is (Ayh)p = — V a Va,/ig + 
{d — l)hp. By considering 

f {V a hP + V r h a )(V a hp + V/3~h a ) > 0, (106) 
Jb 

an integration by parts shows that for all eigenvalues A of Ay, we have A > 2{d — 
1). We have used the fact that the base manifold is Einstein with the same scalar 
curvature as the unit sphere of appropriate dimension. The inequality is saturated 
by modes that are Killing vectors. It is well known that the minimum eigenvalue of 
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Ay on S d is precisely 2(d — 1) [12]. Therefore 



(107) 



The equations for the various modes are found by expanding the linearised equa- 
tions 



-{A L h) ab = c(d+ l)h ab , 



(108) 



into harmonics on the base manifold and considering each linearly independent term 



separately. The various components of (108) may be classified by the number of free 



indices on the base manifold: two, one or none. Each equation can then be written 
so that it is tensorial on the base manifold. There are not many tensors available 
that are linear in the perturbation, they are shown in table 4. 



Table 4: Availab 


e tensors on the base manifold. 


Free base indices 


tensors 


vectors 


scalars 


Two 


h aj3 and (A L h) a p 


V(aV) 


V( a V/3)^ and g a gh 


One 




h a and (Ayh) a 


d a h 


None 






h and Ag/i 



Multiplying these tensors in each equation will be differential expressions involving 
functions of r and t, but the t dependence is just an exponential. The equations 
corresponding to the terms along the diagonal of table 4 will be of the form 



F [drV^drcp 1 , if 1 ] Ah + G [d^tp 1 , drip 1 , p 1 ] h = 0, 



(109) 



where F, G are functionals of the radial functions (p 1 (r). Here / indexes the scalar or 
vector modes because there is more than one mode of each type in the decomposition 
of h a h- In this equation h and A represent a harmonic mode on the base and the 
corresponding differential operator. But we have decomposed into harmonics of these 
differential operators, so Ah = Xh. This then gives a differential equation for the 
p 1 (r) that only depends on the base manifold through the eigenvalue A 



F [d 2 r p T , d^ 1 , p 1 ] X + G [d r V, drp 1 , p 1 ] = 0. 



(110) 



36 



Above the diagonal in table 4 we see that there is only one possible term in each 



case, after noting that V( a Vm/» does not arise from fll08|) . Therefore the correspond- 
ing equations will be of the form 

H[dr<p I ,d r <p I ,<p I ]Vh = Q, (111) 

where Vh is one of the tensors from above the diagonal in table 4. But this simply 
implies the differential equation H = 0, with no dependence on the base manifold. 
Therefore the equations for the perturbations only depend on the base manifold 
through the eigenvalues of the second order differential operators on the base. 

As we saw for the tensor mode, instabilities are associated with the eigenvalue 
being less than some critical value. This is the generic situation and should be 
thought of as being due to a sufficiently negative mass squared. However, in the 
first part of this appendix we saw that the eigenvalues for the scalar and vector 
perturbations were bounded below by the minimum eigenvalues on the sphere of the 
given dimension and unit radius, S . It is expected, but to our knowledge not proven, 
that the standard Schwarzschild-Tangherlini and Schwarzschild-AdS black holes are 
stable against these vector and scalar perturbations and therefore these eigenvalues 
must be larger than any critical value. This suggests that the vector and scalar modes 
do not contribute towards the instability of the generalised Schwarzschild-Tangherlini 
and Schwarzschild-AdS black holes that we have been considering in this work. It 
would be nice, however, to check this explicitly from a perturbation analysis. 



B Transverse tracefree gauge 

He we show that one may impose the tracefree condition in addition to the transverse 
condition for vacuum spacetimes with a cosmological constant. This extends a famil- 
iar argument to the case of nonzero cosmological constant. Suppose the background 
spacetime in d + 2 dimensions satisfies the vacuum field equations, possibly with a 
cosmological constant, 

R ab = c{d + l)g ab , (112) 
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then considering perturbations that are transverse in the sense of ( |l6| ) one may take 
the trace of the perturbed equations 

^(A L h) ab = c(d + l)h ab => nh a a + 2c(d + l)h a a = 0. (113) 

Consider solving the residual gauge freedom equations (|i~7| ) subject to the following 
initial conditions on the spacelike hypersurface, t = to, 

2(V% + V m £m) = -h\, 
2(- V m V m ^o + V m V e m -2c(d+l)£ ) = -V /i a a , (114) 

where m runs over the spatial indices. These equations may always be solved for 
£ and ^ on the hypersurface. This initial data then defines a vector field £ in the 
causal future of the hypersurface using the residual gauge equation (|l7|). We can 
show that this vector field gives a gauge transformation which sets the trace to zero. 
Define f = h a a + 2V a £ a . The initial conditions (|l|) and (0) are seen to imply that 
f = % = at t = to. Now find the equation satisfied by /, using ( |113|) and (|l7|) 

□/ = □/i + 2DV a f a 

= -2c(d + l)h\ + 2V a D£ a - 2c(d + l)V a £ a 

= -2c(d + 1) [/i a a + 2V a £„] = -2c(d + 1)/. (115) 

We see that, perhaps unexpectedly, the equation for / has no source term. This 
fact implies that the vanishing initial conditions force / to vanish identically in the 
causal future of the initial hypersurface. Issues of global hyperbolicity of the initial 
spacetime should not concern us as we are only interested in the future development 
of the perturbation. Thus we have shown that there exists a £ satisfying the residual 
gauge equation ( |17| ) such that the trace of the perturbation may be set to 0. 

C No pure gauge solutions 

We show here that the tensor perturbations considered in this paper cannot be pure 
gauge. Suppose there were a pure gauge mode, h ab = V a £;, + Vft£ a . For the pertur- 
bations we are considering about the metric (|To|), consider the vanishing components 

h n = 2(d r & - I^Ui) = 0. (116) 



38 



Therefore, noting that the form of the perturbation (^) specifies the time depen- 
dence, 

6=5W 1/2 e^6(x). (117) 
Now consider the vanishing component 

froo = 2(d t £o - £tt6) = 0. (118) 
2g(r) 

Together with the form we found for £i this then implies that 

Next consider the vanishing component 

h i = dt£i + drto - ^y£o = 0. (120) 
f{ r ) 

If we now substitute the expressions we have for £o>£i into this equation, we obtain 
a differential equation for g(r) and /(r) that is manifestly not satisfied by any of the 
functions we used in this work. This implies that £o = Ci = 0- Now considering the 
vanishing components 

h 0a = d t i a + d a £ = 0, (121) 

we see that £o = implies that £ a is independent of t. But if £ a is nonzero then it 
must depend on t as e ut in order for the h a p term to have required time dependence. 
The conclusion is thus that £ a = and therefore the perturbations we are considering 
can never be pure gauge. 

In the section on time dependent metrics we took / = and the metric had no 
t components. The argument of the previous paragraph will not work in this case. 
The expression for £i (117) is now just £i = g(r) 1 / 2 ^i(x). We can substitute this into 



the component 



hla = Orta + #a£l " ha = 0. (122) 

r 



Solving for £ a gives 



a 



,S 2 



ds 



dji(x) = K(r)d a ii(x). (123) 
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Finally, substitute these results into the remaining, nonvanishing, component and 
recall the form of the perturbation (^) and the definition of K(r) in ( [L23| ) 

2r 

g{r) 

= 2K(r)V a dp£ 1 (x) + —^g a p(x)l l (x) 

g{r) 1/2 

= r 2 ip{r)K p (x). (124) 

For the last two lines to be equal, we must have either K(r) oc g ^i/ 2 or <x - 
9a0^i of V Q 9/3^i. The former possibility is not true for the functions g(r) that we 
have been considering whilst the latter is not consistent with h a a = 0. Therefore 
there are no pure gauge solutions. We could have used this argument in the previous 
case also, but the argument we gave instead did not use the tracefree property at any 
point. 
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